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Abstract 

There are two sets of four- by- four matrices introduced by Dirac. The first set 
consists of fifteen Majorana matrices derivable from his four 7 matrices. These fifteen 
matrices can also serve as the generators of the group SL(4, r). The second set consists 
of ten generators of the Sp{4) group which Dirac derived from two coupled harmonic 
oscillators. It is shown possible to extend the symmetry of Sp(4) to that of SL(4:, r) 
if the area of the phase space of one of the oscillators is allowed to become smaller 
without a lower limit. While there are no restrictions on the size of phase space 
in classical mechanics, Feynman's rest of the universe makes this Sp(4)-to-SX(4, r) 
transition possible. The ten generators are for the world where quantum mechanics is 
valid. The remaining five generators belong to the rest of the universe. It is noted that 
the groups SX(4, r) and Sp(A) are locally isomorphic to the Lorentz groups 0(3,3) 
and 0(3,2) respectively. This allows us to interpret Feynman's rest of the universe in 
terms of space-time symmetry. 
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1 Introduction 



In 1963, Paul A. M. Dirac published an interesting paper on the coupled harmonic 
oscillators [lj. Using step-up and step-down operators, Dirac was able to construct ten 
operators satisfying a closed set of commutation relations. He then noted that this set 
of commutation relations can also be used as the Lie algebra for the 0(3,2) deSitter 
group applicable to three space and two time dimensions. He noted further that this 
is the same as the Lie algebra for the four-dimensional symplectic group Sp(4). 

His algebra later became the fundamental mathematical language for two-mode 
squeezed states in quantum optics [2j EJ EJ [5]. Thus, Dirac's ten oscillator matrices 
play a fundamental role in modern physics. 

In the Wigner phase-space representation, it is possible to write the Wigner function 
in terms of two position and two momentum variables. It was noted that those ten 
operators of Dirac can be translated into the operators with these four variables [4j [6], 
which then can be written as four-by-four matrices. There are thus ten four-by-four 
matrices. We shall call them Dirac's oscillator matrices. They are indeed the generators 
of the symplectic group Sp(i). 

We are quite familiar with four Dirac matrices for the Dirac equation, namely 
7lj 72 5 73; and 70 • They all become imaginary in the Majorana representation. 

From them we can construct fifteen linearly independent four-by-four matrices. It 
is known that these four-by-four matrices can serve as the generators of the SL(4, r) 
group [61 Ej. It is also known that this SL(4,r) group is locally isomorphic to the 
Lorentz group 0(3,3) applicable to the three space and three time dimensions [Sl[7j. 

There are now two sets of the four- by- four matrices constructed by Dirac. The 
first set consists of his ten oscillator matrices, and there are fifteen 7 matrices coming 
from his Dirac equation. There is thus a difference of five matrices. The question is 
then whether this difference can be explained within the framework of the oscillator 
formalism with tangible physics. 

It was noted that his original 0(3, 2) symmetry can be extended to that of 0(3, 3) 
Lorentz group applicable to the six dimensional space consisting of three space and 
three time dimensions. This requires the inclusion of non-canonical transformations in 
classical mechanics [6]. These noncanonical transformations cannot be interpreted in 
terms of the present form of quantum mechanics. 

On the other hand, we can use this non-canonical effect to illustrate the concept of 
Feynman's rest of the universe. This oscillator system can serve as two different worlds. 
The first oscillator is the world in which we do quantum mechanics, and the second is 
for the rest of the universe. Our failure to observe the second oscillator results in the 
increase in the size of the Wigner phase space thus increasing the entropy [8]. 

Instead of ignoring the second oscillator, it is of interest to see what happens to it. 
It is shown in this paper that Planck's constant does not have a lower limit. This is 
allowed in classical mechanics, but not in quantum mechanics. 

Indeed, Dirac's ten oscillator matrices explain the quantum world for the both 
oscillators. The set of Dirac's fifteen 7 matrices contains his ten oscillator matrices as 
a subset. We discuss in this paper the physics of this difference. 

In Sec. [2j we start with Dirac's four 7 matrices in the Majorana representation 
and construct all fifteen four-by-four matrices applicable to the Majorana form of the 
Dirac spinors. Sec. [3] reproduces Dirac's derivation of the 0(3,2) symmetry with ten 
generators from two coupled oscillators. This group is locally isomorphic to Sp(4), 
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which allows canonical transformations in classical mechanics. 

In Sec. [H we translate Dirac's formalism into the language of the Wigner phase 
space. This allows us to extend the Sp(A) symmetry into the non-canonical region 
in classical mechanics. The resulting symmetry is that of SL(4, r), isomorphic to 
that of the Lorentz group 0(3,3) with fifteen generators. This allows us to establish 
the correspondence between Dirac's Majorana matrices with those SL(A, r) four-by- 
four matrices applicable to the two oscillator system, as well as the fifteen six-by-six 
matrices which serve as the generators of the 0(3,3) group, 

Finally, in Sec. [5l it is shown that the difference between the ten oscillator matrices 
and the fifteen Majorana matrix can serve as an illustrative example of Feynman's rest 
of the universe [8j [9] . 



2 Dirac Matrices in the Majorana Representa- 
tion 

Since all the generators for the two coupled oscillator system can be written as four- 
by-four matrices with imaginary elements, it is convenient to work with Dirac matrices 
in the Majorana representation, where the all the elements are imaginary [TJ [10], [TT] In 
the Majorana representation, the four Dirac 7 matrices are 
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These 7 matrices are transformed like four-vectors under Lorentz transformations. 
From these four matrices, we can construct one pseudo-scalar matrix 

75 = *707i7273 = ( „ ) > ( 2 ) 



-<7 2 



and a pseudo vector 275 7 M consisting of 



(Ti 0\ . .[01 



^7571 =^ a J, ^7572 AJ 

i757o = i(_° J J), 27573 = 2 (7 3 + ° a3 ). (3) 
In addition, we can construct the tensor of the 7 as 

T^ u = - (7^ - 7^) . (4) 
This antisymmetric tensor has six components. They are 

*707i = -2 Q 1 ), 27072 = -2 (" J J), ,7073 = -2 * 3 ) , (5) 



3 



and 



*7i7Si = <(° ^273 = ^( (j 3 " 3 ), ^737i=( CT Q 2 °J. (6) 

There are now fifteen linearly independent four-by-four matrices. They are all 
traceless, their components are imaginary [7]. We shall call these Dirac's Majorana 
matrices. 

In 1963 [T], Dirac constructed another set of four-by-four matrices from two coupled 
harmonic oscillators, within the framework of quantum mechanics. He ended up with 
ten four- by- four matrices. It is of interest to compare his oscillator matrices and and 
his fiftteen Majorana matrices. 

3 Dirac's Coupled Oscillators 

In his 1963 paper [il], Dirac started with the Hamiltonian for two harmonic oscillators. 
It can be written as 

H=~(pl + xl) + 1 ~( P l + xl). (7) 
The ground-state wave function for this Hamiltonian is 

4>o(xi,x 2 ) = -^=exp j-~(a^ + (8) 

We can now consider unitary transformations applicable to the ground-state wave 
function of Eq. ([8]) , and Dirac noted that those unitary transformations are generated 

by m 

Li = - (a\a 2 + a\aij , L 2 = (a\a 2 - a\a\ 

L 3 = ^ ( a i ai ~ a 2 a2 ) ' Ss= 2 
K\ = —— (a\a\ + a\a\ — a^a^ — a 2 a 2 

K 2 = - (a\a{ - a\a\ + a\a\ - a 2 a 2 ^j , 
K i = o ( a i a 2 + ° ia2 ) ' 



a\a\ + a 2 a\ 



Qi = — - (a\a\ — a\ax — a\a\ + 02^2) , 
Q 2 = — - (a\a\ + a\a\ + a\a\ + a 2 a 2 ) , 

Qs = 2 ( a i a 2 ~ ° ia2 ) • ( 9 ) 
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where a) and a are the step-up and step-down operators applicable to harmonic oscilla- 
tor wave functions. These operators satisfy the following set of commutation relations. 



[Li, Lj] — ieijkLk, 
[K i ,K 3 ] = [Q i ,Q j \ 



[Ki,Qj] = -idijSs, 



[Li,Kj] — ieijkKk, [Li,Qj] — itijkQki 



[Ki,S 3 ] 



[Li,S 3 ]=Q, 
-iQi, [Qi,S3,] 



(10) 



Dirac then determined that these commutation relations constitute the Lie algebra 
for the 0(3, 2) deSitter group with ten generators. This deSitter group is the Lorentz 
group applicable to three space coordinate and two time coordinates. Let us use 
the notation (x,y, z,t, s), with (x,y,z) as space coordinates and (t, s) as two time 
coordinates. Then the rotation around the z axis is generated by 



The generators L\ and L2 can be also be constructed. The K% and Q3 will take the 
form 
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(12) 



From these two matrices, the generators K\,K^,Q\,Q2 can be constructed. The gen- 
erator 63 can be written as 



S, 



/o \ 





-i 

\0 i / 



(13) 



The last five-by-five matrix generates rotations in the two-dimensional space of (t,s). 

In his 1963 paper |T|, Dirac states that the Lie algebra of Eq. (fT0|) can serve as the 
four-dimensional symplectic group <Sp(4). In order to see this point, let us go to the 
Wigner phase-space picture of the coupled oscillators. 

3.1 Wigner Phase-space Representation 

For this two-oscillator system, the Wigner function is defined as [H [6] 
W(x 1 ,x 2 ;pi,P2) = (^) J exp{-2i(piyi +P22/2)} 

x ip*(x 1 + yi,x 2 + V2)tp(xi - yi,x 2 - y2)dyidy 2 . (14) 
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Indeed, the Wigner function is denned over the four-dimensional phase space of (xi,pi,X2,P2) 
just as in the case of classical mechanics. The unitary transformations generated by the 
operators of Eq.([9|) are translated into linear canonical transformations of the Wigner 
function [4]. The canonical transformations are generated by the differential opera- 
tors [I]: 



d d \ / d d 

2 \\ 1 dp2 P2 dxi) \ 2 dpi Pl dx 2 
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2{\ 1 dx 2 X2 dxi) \ dp2 P2 dpi 
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2\\ X dp2 P2 dxi) \ 2 dpi Pl dx 2 

Q + l {(x 9 p 9 \ (x 9 p 9 

" ' 2 \ \ dxi dpi J \ dx 2 dp2 

1 (/ a a \ ( a a 

2 LV dpi dxiJ V dp2 dx 2 



1 



d a \ ( a a 



* - --A\ x ^ x ^)-\^^)Y (16) 

3.2 Translation into four-by-four matrices 

For a dynamical system consisting of two pairs of canonical variables x\,pi and X2,P2, 
we can use the coordinate variables defined as [6] 

(r/i,r?2,r/3,?74) = (xi,px, x 2 ,p 2 ) • (17) 

Then the transformation of the variables from iji to is canonical if |12^ H3] 

MJM = J, (18) 

where M is a four-by-four matrix defined by 
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(19) 



According to this form of the J matrix, the area of the phase space for x\ and p\ 
variables remains invariant, and the story is the same for the phase space of X2 and p2- 
we can thren write the generators of the Sp{4) group as 
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These four- by- four matrices satisfy the commutation relations given in Eq. fjlOp . In- 
deed, the deSitter group 0(3,2) is locally isomorphic to the Sp(A) group. The remaining 
question is whether these ten matrices can serve as the fifteen Dirac matrices given in 
Sec. El The answer is clearly No. How can ten matrices describe fifteen matrices? We 
should therefore add five more matrices. 



4 Extension to 0(3,3) Symmetry 

Unlike the case of the Schrodinger picture, it is possible to add five noncanonical 
generators to the above list. They are 
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as well as three additional squeeze operators: 

G --{(x —+x —\ + ( — + — 
1 2 1 \ 1 dx 2 X2 dxi ) \ 1 dp2 P2 dpi 

G -{( — —\ ( — — 

l\\ 1 dp2 P2 dxi) V 2 dpi 1 9x2 

ft = -ilU^+^W^+p^l'- PS) 



2 [\ 9xi dpi J V ^Pl 9x 



These five generators perform well-defined operations on the Wigner function. How- 
ever, the question is whether these additional generators are acceptable in the present 
form of quantum mechanics. 
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In order to answer this question, let us note that the uncertainty principle in the 
phase-space picture of quantum mechanics is stated in terms of the minimum area in 
phase space for a given pair of conjugate variables. The minimum area is determined 
by Planck's constant. Thus we are allowed to expand phase space, but are not allowed 
to contract it. With this point in mind, let us go back to G3 of Eq. ([23|) . which gen- 
erates transformations which simultaneously expand one phase space and contract the 
other. Thus, the G3 generator is not acceptable in quantum mechanics even though it 
generates well-defined mathematical transformations of the Wigner function. 

If the five generators of Eq. (j22[) and Eq. (I23p are added to the ten generators given 
in Eq. (|15p and Eq. (|16p . there are fifteen generators. They satisfy the following set of 
commutation relations. 

[Li, Lj] = ieijkLk, [Si,Sj\ = ieij^Sk, [Li,Sj] = 0, 

[Li,Kj] = ieijkKk, [Li,Qj] = ieijkQk-, [Li-,Gj\ = ieijkGk, 

[Ki,Kj] = [Qi,Qj] = [Qi,Qj] = —ieijkLk, 

[Ki,Qj] = -iSijS 3 , [Qi,Gj] = —iSijSx, [Gi, Kj] = -i5ijS 2 , 

[K h S 3 ] = -iQh [Qh S3] = %K U [G h S 3 ] = 0, 

[K h Si] = 0, [Q h Si] = -iGi, [d, Si] = %Q U 

[K l ,S 2 ] = iG l , [Qi,S 2 }=Q, [Gi,S 2 } = -iKi. (24) 



As we shall see in Subsec. 14.21 this set of commutation relations serves as the Lie 
algebra for the group SL(4, r) and also for the 0(3,3) Lorentz group. 

These fifteen four-by-four matrices are written in terms of Dirac's fifteen Majorana 
matrices, and are tabulated in Table [TJ There are six anti-symmetric and nine sym- 
metric matrices. These anti-symmetric matrices were divided into two sets of three 
rotation generators in the four-dimensional phase space. The nine symmetric matrices 
can be divided into three set of three squeeze generators. However, this classication 
scheme is easier to understand in terms the group 0(3,3), discussed in Subsec. [4^21 

4.1 Non-canonical Transformations in Classical Mechan- 
ics 

In addition to Dirac's ten oscillator matrices, we can consider the matrix 

*=KJ -J- (25) 

which will generate a radial expansion of the phase space of the first oscillator, while 
contracting that of the second phase space p3], as illustrated in Fig. [TJ What is 
the physical significance of this operation? The expansion of phase space leads to an 
increase in uncertainty and entropy [HI Q3] . 

The contraction of the second phase space has a lower limit in quantum mechanics, 
namely it cannot become smaller than Planck's constant. However, there is not such 
lower limit in classical mechanics. We shall go back to this question in Sec. [5j 

In the meantime, let us study what happens when the matrix G3 is introduced into 
the set of matrices given in Eq. ([20j) and Eq. ([2Tj) . It commutes with S 3 ,L 3 ,K\,K 2 ,Q\, 



Table 1: SL(A, r) and Dirac matrices. Two sets of rotation generators and three sets of boost 
generators. There are 15 generators. 
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Third component 


Rotation 
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and Q2- However, its commutators with the rest of the matrices produce four more 
generators: 

[G 3 ,L 1 ] = iG 2 , [G 3 ,L 2 ] = -iG 1 , [G 3 ,K 3 ]=iS 2 , [G 3 ,Q 3 ] = -iS u (26) 
where 

i [0 I\ 1 / -a 2 

If we take into account the above five generators in addition to the ten generators of 
Sp(4), there are fifteen generators. These generators satisfy the set of commutation 
relations given in Eq. (|24j) . 

Indeed, the ten Sp(4) generators together with the five new generators form the 
Lie algebra for the group SL(A, r). There are thus fifteeen four-by-four matrices. They 
can be written in terms of the fifteen Majorana matrices, as given in Table [U 

4.2 Local Isomorphism between 0(3,3) and SL(4,r) 

It is now possible to write fifteen six-by-six matrices which generate Lorentz transfor- 
mations on the three space coordinates and three time coordinates [6]. However, those 
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Figure 1: Expanding and contracting phase spaces. Canonical transformations leave the 
area of each phase space invariant. Non-canonical transformations can change them. Yet 
the product of these two areas remain invariant. 



Table 2: Three- by-three matrices constituting the two- by-two representation of generators 
of the 0(3, 3) group. 
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matrices are difficult to handle and do not show existing regularities. In this section, 
we write those matrices as two- by- two matrices of three- by-three matrices. 

For this purpose, we construct four sets of three- by-three matrices given in Table [2j 
There are two sets of rotation generators 

applicable to the space and time coordinates respectively. 

There are also three sets of boost generators. In the two-by-two representation of 
the matrices given in Table [21 they are 

where the three- by-three matrices Ai,Bi,Ci, and Di are given in TableEJ and Ai,Bi,Ci, 
are their transposes respectively. 

There is a four-by-four Majorana matrix corresponding to each of hese fifteen six- 
by-six matrices, as given in Table [TJ 

There are of course many interesting subgroups. The most interesting case is the 
0(3, 2) subgroup, and there are three of them. Another interesting feature in that there 
are three time dimensions. Thus, there are also 0(2,3) subgroups applicable to two 
space and three time coordinates. This symmetry between space and time coordinates 
could be an interesting future investigation. 



5 Feynman's Rest of the Universe 

In his book on statistical mechanics [9], Feynman makes the following statement. When 
we solve a quantum-mechanical problem, what we really do is divide the universe into 
two parts - the system in which we are interested and the rest of the universe. We then 
usually act as if the system in which we are interested comprised the entire universe. 
To motivate the use of density matrices, let us see what happens when we include the 
part of the universe outside the system. 

We can use two coupled harmonic oscillators to illustrate what Feynman says about 
his rest of the universe. One of the oscillators can be used for the world in which we 
make physical measurements, while the other belongs to the rest of the universe |8j. 

Let us start with a single oscillator in its ground state. In quantum mechanics, 
there are many kinds of excitations of the oscillator, and three of them are familiar to 
us. First, it can be excited to a state with a definite energy eigenvalue. We obtain 
the excited-state wave functions by solving the eigenvalue problem for the Schrodinger 
equation, and this procedure is well known. 

Second, the oscillator can go through coherent excitations. The ground-state oscil- 
lator can be excited to a coherent or squeezed state. During this process, the minimum 
uncertainty of the ground state is preserved. The coherent or squeezed state is not 
in an energy eigenstate. This kind of excited state plays a central role in coherent 
and squeezed states of light which have recently become a standard item in quantum 
mechanics. 

Third, the oscillator can go through thermal excitations. This is not a quantum 
excitation, but is a statistical ensemble. We cannot express a thermally excited state 
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by making linear combinations of wave functions. We should treat this as a canonical 
ensemble. In order to deal with this thermal state, we need a density matrix. 

For the thermally excited single-oscillator state, the density matrix takes the form [9j 

mils]. 

~ (30) 



p{x,y) = ( 



1 



-1/T 



where the absolute temperature T is measured in the scale of Boltzmann's constant, 
and <pk{x) is the k-th excited state wave oscillator wave function. The index ranges 
from to oo. 

We also use Wigner functions to deal with statistical problems in quantum mechan- 
ics. The Wigner function for this thermally excited state is [9j HI [15] 



W T (x,p) = - [ e 

TT J 



■2ipz 



which becomes 



W T 



tanh(l/2T) 



TT 



exp 



This Wigner function becomes 



W 



TT 



■ cxp 



p(x — z, x + z)dz, 



(x 2 +p 2 ) tanh(l/2T) 



(x 2 + p 2 



(31) 



(32) 



(33) 



when T = 0. As the temperature increases, the radius of this Gaussian form increases 
from one to [14| 

1 (34) 
Vtanh(l/2T) 

The question is whether we can derive this expanding Wigner function from the 
concept of Feynman's rest of the universe. In their 1999 paper [8], Han et al. used 
two coupled harmonic oscillators to illustrate what Feynman said about his rest of 
the universe. One of their two oscillators is for the world in which we do quantum 
mechanics and the other is for the rest of the universe. However, these authors did not 
use canonical transformations. In Subsec. 15. 1\ we summarize the main point of their 
paper using the language of canonical transformations developed in the present paper. 

Their work was motivated by the papers by Yurke et al. |17j and by Ekert et al. |18j , 
and the Barnett-Phoenix version of information theory [19]. These authors asked the 
question of what happens when one of the photons is not observed in the two-mode 
squeezed state. 

In Subsec. 15. 2\ we introduce another form of Feynman's rest of the universe, based 
on non-canonical transformations discussed in the present paper. For a two-oscillator 
system, we can define a single-oscillator Wigner function for each oscillator. Then non- 
canonical transformations allow one Wigner function to expand while forcing the other 
to shrink. The shrinking Wigner function has a lower limit in quantum mechanics, while 
there is none in classical mechanics. Thus, Feynman's rest of the universe consists of 
classical mechanics where Planck's constant has no lower limit. 

In Subsec. 15.31 we translate the mathematics of the expanding Wigner function into 
the physical language of entropy. 
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Figure 2: Two-dimensional Gaussian form for two-coupled oscillators. One of the variables 
is observable while the second variable is not observed. It belongs to Feynman's rest of the 
universe. 



5.1 Canonical Approach 

Let us start with the ground-state wave function for the uncoupled system. Its Hamil- 
tonian is given in Eq.([7]), and its wave function is 



i>o(x 1> x 2 ) 



1 



7T 



exp 



2 ( X? + 



(35) 



We can couple these two oscillators by making the following canonical transformations. 
First, let us rotate the coordinate system by 45° to get 



1 / \ 1 , 

-j={Xi+X 2 ), -j=(Xi-X 2 ), 

Let us then squeeze the coordinate system: 

-j= [xx + x 2 ) , -j=[xx-x 2 ). 
Likewise, we can transform the momentum coordinates to 



(36) 



(37) 



-j=(j>l+P2), -j=[px-P2)- 



(38) 



Equations (|37p and (|38p constitute a very familiar canonical transformation. The 
resulting wave function for this coupled system becomes 



ip v (xx,x 2 ) = -j=exp \-\ \e v (xx - x 2 ) 2 + e v (xx + x 2 ) 2 
This transformed wave function is illustrated in Fig. [2l 



(39) 
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As was discussed in the literature for several different purposes [H [2IJ [22j [23] , this 
wave function can be expanded as 

1 / r]\ k 

^ v (xi,x 2 ) = —^^2[tanh-J fa^fafa), (40) 

where the wave function 4>k4>(x) and the range of summation are defined in Eq. (l30p . 
From this wave function, we can construct the pure-state density matrix 

p(x 1 ,x 2 ;x' 1 ,x' 2 ) = Tp, n (x 1 ,x 2 )?p T1 (x' 1 ,x 2 ), (41) 

which satisfies the condition p 2 = p: 

p(xi,x 2 ;x' 1 ,x' 2 ) = / p(xi,x 2 ]x'{,X2)p(x'{,x 2 ]x' 1 ,x' 2 )dx'{dx 2 . (42) 



If we are not able to make observations on the x 2 , we should take the trace of the 
p matrix with respect to the x 2 variable. Then the resulting density matrix is 

p(x,x') = J ifJ v (x,x 2 ) {if) v (x',x 2 )}* dx 2 . (43) 

Here, we have replaced x\ and x\ by x and x' respectively. If we complete the integra- 
tion over the x 2 variable, 



p(x, x') = ( l —— \ exp 

\7rcosn77 / 



(x + x') + (x — x') cosh rj 



4 cosh r] 

The diagonal elements of the above density matrix are 



(44) 



p(x,x) = ( — ^ exp (—x 2 / cosh rj). (45) 

V 71" cosh rj J ^ ' 

With this expression, we can confirm the property of the density matrix: Tr(p) = 1. 
As for the trace of p 2 , we can perform the integration 

Tr(p 2 \= ( p(x, x')p(x', x)dx'dx = , (46) 

^ ' J cosh 77 

which is less than one for nonzero values of r\. 

The density matrix can also be calculated from the expansion of the wave function 
given in Eq. (|40p . If we perform the integral of Eq. (|43p . the result is 



P( x > x ') = ( cosh^/2) ) ^ ( tanh l) ^{xWkix'), (47) 

which leads to Tr(p) = 1. It is also straightforward to compute the integral for Tr(p 2 ). 
The calculation leads to 

1 \* / n\ 4 * 



Tr M = Ub)J ?( tanh f) ■ <48) 

The sum of this series becomes to (1/ cosh r/), as given in Eq. (|46p . 
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We can approach this problem using the Wigner function. The Wigner function for 
the two oscillator system is [4] 



W (xi,pi;x 2 ,P2) 



— exp 

7T 



x\ + p\ + A + pl 



(49) 



If we pretend not to make measurement on the second oscillator coordinate, the X2 and 
P2 variables have to be integrated out [8]. The net result becomes the Wigner function 
for the first oscillator. 

The canonical transformation of Eq. (|37p and Eq. (|38|) changes this Wigner function 

to 



W(xi,x 2 ;pi,P2) 



cxp 



e n {xi -x 2 f + e~ 11 (x 1 + x 2 f 



+e~ r '( Pl - p 2 ) 2 + e"(pi + p 2 ) 2 ] } • 



(50) 



If we do not observe the second pair of variables, we have to integrate integrate this 
function over %2 and p2- 



(51) 



W v (xi,p 1 ) = J W(x 1 ,X2;pi,P2)dx 2 dp2, 

and the evaluation of this integration leads to [8] 

1 



W v (x,p) 



IT COsh 7] 



exp 



x 2 + p 2 

COsh 7] 



(52) 



where we use x and p for x\ and p\ respectively. 

This Wigner function is of the form given in Eq. (|32p for the thermal excitation, if 
we identify the squeeze parameter rj as [20J 



cosh 7] 



1 



tanh(l/2T)' 



(53) 



The failure to make measurement on the second oscillator leads to the radial expansion 
of the Wigner phase space as in the case of the thermal excitation. 



5.2 Non-canonical Approach 

As we noted before, among the fifteen Dirac matrices, ten of them can be used for 
canonical transformations in classical mechanics, and thus in quantum mechanics. 
They play a special role in quantum optics [21 [31 [H [5] . 

The remaining five of them can have their roles if the change in the phase space 
area is allowed. In quantum mechanics, the area can be increased, but it has a lower 
limit called Plank's constant. In classical mechanics, this constraint does not exist. 
The mathematical formalism given in this paper allows us to study this aspect of the 
system of coupled oscillators. 

Let us choose the following three matrices from those in Eq. (|20p and Eq. (|2ip . 
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They satisfy the closed set of commutation relations: 



[S 3 ,K 2 ]=iQ 2 , [S 3 ,Q 2 ] = -iQ 3 , [K 2 ,Q 2 ] = -iS 3 . (55) 

This is the Lie algebra for the Sp(2) group, This is the symmetry group applicable 
to the single-oscillator phase space [3J, with one rotation and two squeezes. These 
matrices generate the same transformation for the first and second oscillators. 
We can choose three other sets with similar properties. They are 



and 



These matrices also satisfy the commutation relations given in Eq. (|55p . In this 
case, the squeeze transformations take opposite directions in the second phase space. 

Since all these transformations are canonical, they leave the area of each phase 
space invariant. However, let us look at the non-canonical generator G 3 of Eq. (|25p . It 
generates the transformation matrix of the form 



e 7 * 
e-* 



(59) 



If r\ is positive, this matrix expands the first phase space while contracting the second. 
This contraction of the second phase space is allowed in classical mechanics, but it has 
a lower limit in quantum mechanics. 

The expansion of the first phase space is exactly like the thermal expansion result- 
ing from our failure to observe the second oscillator which belongs to the rest of the 
universe. If we expand the system of Dirac's ten oscillator matrices to the world of 
his fifteen Majorana matrices, we can expand and contract the first and second phase 
spaces without mixing them up. We can thus construct a model where the observed 
world and the rest of the universe remain separated. In the observable world, quantum 
mechanics remains valid with thermal excitations. In the rest of the universe, since the 
area of the phase space can become small without lower limit, only classical mechanics 
is valid. 

During the expansion/contraction process, the product of the areas of the two phase 
space remains constant. This may or may not be an extended interpretation of the the 
uncertainty principle, but we choose not to speculate further on this issue. 

Let us turn our attention to the fact that the groups SL(4, r) and Sp(4) are locally 
isomorphic to 0(3, 3) and 0(3,2) respectively. This means that we can do quantum 
mechanics in one of the 0(3,2) subgroups of 0(3, 3), as Dirac noted in his 1963 pa- 
per [T]. The remaining generators belong to Feynman's rest of the universe. 



5.3 Entropy and the Expanding Wigner Phase Space 

We have seen how Feynman's rest of the universe increases the radius of the Wigner 
function. It is important to note that the entropy of the system also increases. 
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Let us go back to the density matrix. The standard way to measure this ignorance 
is to calculate the entropy defined as [TBI IM l25l l26l 127] 



S = -Tr(phx(p)), (60) 

where S is measured in units of Boltzmann's constant. If we use the density matrix 
given in Eq. (|44j) . the entropy becomes 

S = 2 jcosh 2 H^j In (^cosh |^ - sinh 2 f|J In fsinh ^ | . (61) 

In order to express this equation in terms of the temperature variable T, we write 
Eq.([53]) as 

1 + e~ x l T 

cosh 77 = - r , (62) 



.2 M 



1 . o /n\ e" 1 / 7 



which leads to 

crf If) = TTFW- srf (I) = TTFW' (63) 

Then the entropy of Eq. (|6ip takes the form [S] 

S=(^) i * A - ^ fl - e" 1 ^ • (64) 



Ty lexp(l/T) - 1 

This familiar expression is for the entropy of an oscillator state in thermal equilibrium. 
Thus, for this oscillator system, we can relate our ignorance of the Feynamn's rest of 
the universe, measured by of the coupling parameter ry, to the temperature. 



6 Concluding Remarks 

In this paper, we started with the fifteen four-by-four matrices for the Majorana repre- 
sentation of the Dirac matrices, and the ten generators of the Sp(4) group correspond- 
ing to Dirac's oscillator matrices. Their explicit forms are given in the literature [3 [7], 
and their roles in modern physics are well known [31 HI [IT] . We re-organized them into 
tables. 

The difference between these two representations consists of five matrices. The 
physics of this difference is discussed in terms of Feynman's rest of the universe [9]. 
According to Feynman, this universe consists of the world in which we do quantum 
mechanics, and the rest of the universe. In the rest of the universe, our physical laws 
may or may not be respected. In the case of coupled oscillators, without the lower 
limit on Planck's constant, we can do classical mechanics but not quantum mechanics 
in the rest of the universe. 

In 1971, Feynman et al. [30] published a paper on the oscillator model of hadrons, 
where the proton consists of three quarks linked up by oscillator springs. In order to 
treat this problem, they use a three-particle symmetry group formulated by Dirac in 
his book on quantum mechanics |31} [32] . An interesting problem could be to see what 
happens to the two quarks when one of them is not observed. Another interesetinq 
question could be to see what happens to one of the quarks when two of them are not 
observed. 
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Finally, we note here that group theory is a very powerful tool in approaching 
problems in modern physics. Different groups can share the same set of commutation 
relations for their generators. Recently, the group SL(2, c), through its correspondence 
with the <SO(3, 1) has been shown to be the underlying language for classical and 
modern optics [UGH]. In this paper, we exploited the correspondence between SX(4, r) 
and 0(3,3), as well as the correspondence between Sp(4) and 0(3,2) which was first 
noted by Paul A. M. Dirac pQ. 

There could be more applications of group isomorphisms in the future. A compre- 
hensive list of those correspondences is given in Gilmore's book on Lie groups [29]. 
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